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Introduction - Atom interferometers offer an 

unprecedented precision in inertial measurements. Sup- 
plemented with a highly coherent input source provided 
by Bose-condensed atoms 0, la H H UH > atom inter- 
ferometers may potentially supersede the conventional 
laser-based devices. In a recent experiment |ll| a minia- 
ture Michelson- type interferometer was realized on an 
atom chip HI El El 13 13 113: thus further approach- 
ing practical implementations of the device. Generally, 
miniaturization of atomic devices leads to an increased 
role of interatomic interactions, due to higher densities 
and density gradients. Indeed, a strong suppression of 
contrast was observed in E3 f° r longer durations of the 
interferometric cycle. The goal of our paper is to ex- 
plain this effect and suggest recipes for suppressing the 
interaction-related fringe degradation. 

The role of interatomic interactions in interferometric 
processes has been studied by several authors 0, 0, |3 
l2lj , with the main emphasis on the potential loss of first- 
order coherence. In our paper we focus on a different 
effect: distortion of the interferometric path due to the 
mean-field pressure. 

Interferometric scheme.- In this paper we consider the 
Michelson interferometric scheme (see Fig. ma)). Atoms 
are supposed to be confined transversally by a monomode 
atom guide, with no transverse excitations allowed. The 
initial state of atoms is a perfectly coherent state 
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normalized for convenience to the total number of atoms 
N: J_™dz\ X (z)\ 2 = N. We further assume that at 
every stage of the process the wave function can be ap- 
proximately decomposed into a sum of three spatial har- 
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ing the interferometric process, we show below that they 
are (a) small under typical conditions, and (b) if neces- 
sary can be taken into account a posteriori.) Splitting, 
reflection, and recombining pulses perform the following 
instant transformations of the vector (<E>_i, <E>o, 
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Such ideal interferometric elements were proposed in |2^ | 
and successfully experimentally realized in II]. The 
splitting, reflection, and recombination pulses are applied 
in succession, separated by an equal time interval T. Im- 
mediately after the recombination pulse the population 
of atoms in the central peak is detected; this constitutes 
the interferometric signal. 

Between the splitting and recombination pulses the 
wave function can be approximately decomposed as 
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where the phases 
real. Here 



will be shown to be approximately 
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for < t < T 



±V Q T =f Vq (t — T) for T <t< 2T 
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are the classical trajectories, originating at z = 0, corre- 
sponding to the right (+) and left (-) arms of the inter- 
ferometer; 



where $ n (z, t) are slow functions of coordinate. (Note 
that even though higher harmonics can be generated dur- 
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FIG. 1: Interferometric loop of the Michelson interferometer 
with single atoms (a) and with Bose condensates (b). 



are the corresponding classical velocities; finally, e+(i) = 
e_(i) = Eq = mVq/2 are the corresponding kinetic en- 
ergies. The velocity Vq is given by Vq = hQ/m, where 
to is the atomic mass. 

The resulting differential phase shift 



A0(z) = 0+(z, 2T)-cf>-(z, 2T) 



(7) 



consists of two parts: A<f>(z) = A0 signal + A0 distortion (z). 
The first, spatially independent part is the useful sig- 
nal, related to the effect the interferometer measures. 
The second part is the result of the distortion caused by 
unaccounted-for external fields and, in our case, mean- 
field interactions. The distortion phase shift leads to 
two effects. The first is a correction to the signal phase 
shift. This effect can in principle be accounted for, if 
the nature of the distortion is known. The second ef- 
fect is a degradation of contrast in the interferometric 
signal. This degradation can not be eliminated easily if 
A0distortion(^) changes substantially over the length of 
the atomic cloud. If both effects are taken into account, 
then the interferometric signal 

5(A0 signa i) = N- 1 x 

Xj 
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can be shown to be 

1 M 

^(A^signai) = - + — COs(A0 signa i - 5) , 



where the fringe contrast M and fringe shift 5 are given 
by 



M = yjA 2 + B 2 
5 = arg(A B) , 
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and 
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The goal of this paper is to calculate the fringe contrast 
degradation caused by mean-field interactions, compare 
the results with the experimentally observed values, and 
suggest methods for eliminating this degradation. 

Analysis of the mean-field effects- We describe the 
evolution of the atomic cloud using the nonlinear 
Schrodinger equation 



2m dz' 



-4j = U(z, t)1> + sidM V 
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where U(z, t) is an external field comprising the ob- 
ject of the measurement and any other auxiliary fields 
present; = —?i 2 //i am is the one-dimensional cou- 
pling constant; am = (— a 2 L /2a) [1 — C(a/a±_)] is the one- 



dimensional scattering length 
mass; a± 



/i = to/2 is the reduced 



y/h/(UjJ± is the size of the transverse ground 
state of the guide; C = 1.4603 . . . (see |23); lastly, a is the 
three-dimensional s-wave scattering length. We further 
decompose the wave function into a quasi- Fourier series 
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where each term obeys 
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dt 2m dz 2 



The initial condition for the system l|12|l is given by the 
result of applying the beamsplitter @ to the initial con- 
dition QJ. We get 

^V±i(M = 0+) = . 

Let us assume that the i))±\ components of the wave 
function remain dominant, as they are initially, and then 
verify this assumption for self-consistency. Here and 
throughout the paper we will suppose that the beamsplit- 
ter recoil energy Eq is much greater than the mean-field 
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0) can be traced to the difference between the Hartree- 
type interaction among atoms in the condensate, and the 
Hartree-Fock-type interaction between atoms in the con- 
densate and those out of the condensate. 

In what follows we study the effect of the mean-field 
interaction l|16|) on the differential phase shift (7J and 
fringe contrast ©• 

Differential phase shift - From now on we will be using 
the representation |@J for the wave function, where the 
phase factors 4>±(z, t) (generally complex) are still un- 
known. The equation (|15|l under the representation Q 
leads to the following equations for the phase factors 



d 

— <j)± = -U±(z, t)/h + D±(z, t) 



(17) 



FIG. 2: An artistic view of the mean-field effects in a Michel- 
son interferometer. Notice that the mean-field potential is 
different for the right arm (a) and the left arm (b). 



Under these assumptions the strongest higher harmonics 
generated are 



Here 



^±3 = - 
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and it is indeed small provided the parameter (|13[) is 
small. Another assumption used to derive (|14f> was that 
the beamsplitter recoil energy Eq dominates another en- 
ergy scale associated with the deviation of the momen- 
tum distribution of the dominant harmonics i])±x from 
that of strict 5-peaks around p = ±h~Q. As we will see 
from the following, in a realistic experiment this energy 
scale satisfies an even stricter requirement of being dom- 
inated by the "time-of-flight" energy h/T. 

In general, according to our findings the dynamics 
of higher harmonics is simply reduced to an adiabatic 
following of the dominant if)±i ones. This observation 
greatly simplifies further analysis. In short, one may 
safely exclude the couplings to higher harmonics from 
the equations of evolution for ip±\, obtain results, and 
add the correction (|14l) a posteriori. 

The equations for the dominant harmonics ip±i now 
become 



where the effective potentials for the right (+) and left 
(-) arms of the interferometer read 

U ± (z, t) = U(z, t) + ^\i> ±1 \ 2 + g 1D \i> Tl \ 2 . (16) 

The first part of the mean-field potential is the contribu- 
tion of the part of the cloud to which the given interfer- 
ometer arm belongs. The second, twice-as-strong part is 
the influence of the opposite arm. The factor of two dif- 
ference between the two contributions (illustrated in Fig. 
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the classical velocity v± is given in JSJl; the potentials 
((To) are 



U±(z, t) = U(z, t) 



(19) 



the classical coordinate z T (t) is represented by Fi- 
nally, the correction term D is 



D± = D.; : ± + -Ddyn., ± + -Ddyn.-init., ± 

where 



D 



init.. ± 



D 



j,± 



D, 



2m dz 
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(21) 

(22) 
(23) 



and 



D 



init., ± 



1 



X(z-z±{t)) 



2m x( z - z ±{t)) 
Now we decompose the phase shift 

4>± = 4>± + 5<t>± 

into a sum of the principal part obeying 

JL$ ± = -u(z,t)±/h , 



(24) 



(25) 



(26) 



and the correction originating from the D term in \17f. 
We assume that the correction 5(f>± is small, 



54>± < 1 
which we will justify later. 



(27) 
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Finally the differential phase shift reads 

Acj)(z) = -fT 1 / dt' x (28) 
Jo 

{U + (z-z+(2T) + z + (t'),t') 
-U-(z-Z-(2T)+z-(t'), t')} 

In what follows we will assume that the initial atomic 
wave function has the Thomas-Fermi profile 

X (z) = ^f^(R 2 -z 2 )9(R-\z\) (29) 

Having in mind certain applications we will also add an 
external harmonic oscillator potential 

U(z, t) = ^ (30) 

of an arbitrary frequency loq. We will further classify 
both as distortion effects and set the useful signal to zero. 
To the order of approximations we made, the caused by 
distortion fringe shift and fringe contrast degradation an- 
alyzed below are independent of the useful signal, and 
thus the latter can indeed be omitted. The distortion 
differential phase shift now becomes equal to the total 
one: A0 dist ortion(z) = Acj>(z). 

Interaction-induced loss of contrast: small 
interferometers- First consider the case of small 
interferometers, where the interferometric arms are sub- 
stantially overlapped during the whole interferometric 
cycle, VqT <C R. In this case we can easily compute the 
combined mean-field and harmonic-oscillator contribu- 
tion to the differential phase shift Ac/> (see , as well 
as the resulting fringe contrast M (see [HJl . They read 

A</> d istortion(z) = 2Kz (31) 

and 

M=\f(2KR)\ . (32) 

Here 

K = mCj 2 T 2 V Q /h (33) 

has the meaning of the half of the differential momentum 
acquired by atoms in the mean field and in the harmonic 
oscillator field (see Fig. ^b)), 

£ 2 = 2^ BEC -^ , (34) 

and 

/(£)=3(sin(0-£cos(£))/£ 3 . (35) 

Here and below c^bec is the frequency of the harmonic 
trap for which the state (|29|l would be the ground state, 
and R is the Thomas-Fermi radius. Note that for the 
configuration considered the fringes are suppressed, but 
either not shifted at all or shifted by 7r (see (0). Notice 



also that for a particular case of ujq = a/^bec the dif- 
ferential phase shift vanishes and thus the fringe contrast 
is strictly 100%. We will discuss the implications of this 
phenomenon below. 

The expression (|32|) for the fringe contrast in small 
interferometers is the first principal result of this paper. 

Remedies for the contrast degradation: small 
interferometers- As one can see from (|31H . the 
distortion of the differential phase shift in small inter- 
ferometers disappears completely if the frequency of the 
external harmonic trap is chosen to be \[2 higher than 
the frequency of the mean-field potential: 

^0 = \/2w B EC =>■ A0 dist ortion = . (36) 

This can be realized in two ways: 

(a) The first scheme is applicable if the longitudinal fre- 
quency uj can be controlled independently of the trans- 
verse frequency uj±. In this case, one should start by 
preparing the condensate in the ground state of the lon- 
gitudinal trap. Then, just before the splitting pulse, one 
should increase the longitudinal frequency by a factor of 

in a short ramp (Fig. |^a)). 

(b) The second scheme assumes that both longitudi- 
nal and transverse frequencies are controlled by the same 
source, and thus the ratio between them is always the 
same. Then the change in the longitudinal frequency 
will affect the nonlinear coupling constant g (to which the 
square of the condensate frequency w B ec i s linearly pro- 
portional) due to the simultaneous change in the trans- 
verse frequency. Note that g is linearly proportional to 
lu±. In this case one satisfies the condition l|3fc>[) by in- 
creasing, prior to the splitting pulse, both the longitudi- 
nal and the transverse frequencies by a factor of 2 (Fig. 
lib)). 

Note that both schemes are completely insensitive to 
possible fluctuations in the number of particles in the 
condensate. 

Interaction-induced loss of contrast: large 
interferometers- Consider now the the case of large 
interferometers where for some period of time during the 
cycle the arms are totally spatially separated: VqT > R. 
The differential phase shift in this case reads 

Acj) distmtion (z) = 2K z- z 3 /l 3 , (37) 

where 

K = -mulT 2 V Q /h (38) 

/ = (7710^0/3^)* . (39) 

Notice that unlike in the small interferometer case, for 
no choice of parameters the differential phase shift can 
be completely eliminated, and that the contribution from 
the harmonic potential grows with the time duration of 
the cycle while the one from the mean-field is stationary. 
This makes us to believe that for large interferometers 
the most promising implementation will be the free-space 
one with no longitudinal confinement present. 
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FIG. 3: Two schemes for the preparation of the initial wave 
packet, designed to eliminate the interatomic-interaction- 
induced degradation of the fringe contrast in small interfer- 
ometers, (a) Situation when the longitudinal confinement can 
be controlled independently from the transverse one. (b) Sit- 
uation when the longitudinal and transverse confinements are 
linearly linked one to another. 




FIG. 4: Fringe contrast ratio vs. the universal parameter r\ 
for large interferometers. 



In the case of ojq = the fringe contrast assumes the 
following compact expression: 

M = F{ti) , (40) 

(see Fig. where 

77 = (2rmjl EC R 3 /3hV Q )i 

= (gi D N/hV Q )i , (41) 
is the parameter governing the fringe contrast, 

T{r,) = ~iF 2 (l/6; 1/2, 7/6; ~n 6 /16) 

- sin(77 3 /2)/?7 3 , (42) 

and n F m (ai, . . . , a n ; bi, . . . , & m ; £) is the generalized 
hypergeometric function. 

Notice that the parameter rj depends on neither the 
shape of the atomic cloud nor on the interferometric cy- 
cle duration. For a given set of atomic and waveguide 
parameters the large contrast requirement 77 <C 1 defines 
a universal limit for the number of atoms: 

N<nV Q / 9lIi , (43) 

The expression (|4U|) is the second principal result of our 
paper. 

Limits of the validity of our computational scheme- 
In order for our conclusions be valid the correction 127|) 
originating from the neglected terms in the kinetic energy 
must be small. We have performed a thorough investi- 
gation aimed at understanding the physical meaning of 
the neglected corrections and estimating their value. The 
results are as follows. 

The correction 5(/)± can be decomposed into a sum of 
four terms 

8<f>± = S(j).l ± + <5</>dyn.. ± + 

<50dyn.-init., ± + <50init., ± (44) 



with the following interpretation: 

e ^,,± K \j u± \-h (45) 

g<50dyn., ±+<50dyn.-init., ±+50init., ± (^g^J 

^{X{z~z±{t))\e-^> t(pu ^ )/2m %{z--z ± m- 

The first correction is related to the expansion factor ( Ja- 
cobian) of the bunch of trajectories of classical particles 
moving in the field U±. The second, third, and forth 
corrections originate from the neglected kinetic energy, 
both the initial kinetic energy (coming from the momen- 
tum distribution of x) and that acquired in the field U± . 
The above corrections (together with the nonlinear cor- 
rection H13|) ~) lead to the following requirements for the 
validity of the approximations used: 

ej = (CoT) 2 « 1 (47) 

e dyn . = mu A z 2 T 3 /h < 1 (48) 

e init = hT/mR 2 < 1 (49) 

e„i = (u>R/Vq) 2 « 1 , (50) 

where z is the typical atomic coordinate, given by z ~ 
R (z ~ VqT) for small (large) interferometers. For a 
typical set of parameters of the JILA experiment with 
u> = 2n X 3.2Hz and T — 10~ 3 s (corresponding to the 
small interferometer case) , the values of these parameters 
are indeed small, validating our approximation: Ej = 

£dyn.-init. = 4.1 X 10~ 4 , Edyn. = 4.7 X 10~ 4 , £ in it. = 3.6 X 

10" 4 , and £ n i = 3.6 x 10" 3 . 

Comparison with the JILA experiment- The parame- 
ters of the JILA's experiment on Michelson interferome- 
ter on a magnetic chip [ll| lie in the range intermediate 
between the small and large interferometer regimes and 
requires no assumption on the distance between the arms 
Vq vis a vis the cloud size R. The equation for the fringe 
contrast in this case must be integrated numerically. In 
[ll| a time- and space-localized pulse of magnetic field 
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was used as the phase signal. A stationary harmonic po- 
tential of frequency uj = 2t: x 5 Hz was present in each 
realization. The contrast was traced as a function of 
the duration of the interferometric cycle 2T, as depicted 
in Fig. [3] Other parameters read ui± = 2tt x 100 Hz, 
Q = 4ir/A, where A = 780 nm is the wavelength of light 
used to produce the interferometric elements, R = 45 /im, 
and a = 100.4 q,b- In the experiment the number of atoms 
varied from one value of the cycle duration to another; 
these numbers are shown in Fig. |Sf a) . The value of wbec 
was extracted from w BEC = (4/9) (?id A/mi? 3 . The re- 
sults of the comparison are shown at the Fig. [Sfa). 

One can see from the Fig. that for the parameters 
chosen the role of interatomic interactions in contrast 
degradation is relatively small and the main source of 
the effect is the stationary harmonic trap. This is en- 
tirely unexpected since the strength of the interactions 
was very close to the strength of the trap, typically 
w bec = -25-^.4u;g, and moreover in the small interferom- 
eter regime the strength interactions become multiplied 
by a factor of 2 of Fock origin (see (|34|) ■ One can show 
further that the relatively weak role of interactions in the 
JILA experiment is not related to any small parameter, 
but is solely an interplay of numerical prefactors. To il- 
lustrate this point we show at the Fig. [3b) a theoretical 
prediction for N = 4.5 x 10 4 atoms exhibiting a dominant 
role of the interatomic interactions. 

Summary and outlook — 

(1) We have developed a simple computational scheme 
that allows to include the effect of interatomic interac- 
tions in calculation of fringe shift and fringe degradation 
in waveguide based atom interferometers. 

(2) In two cases we have found simple analytic expres- 
sion for fringe contrast. These cases are: small interfer- 
ometers where the spatial separation between the arms 
is much smaller than the atomic cloud size, VqT <C R; 
and large interferometers where at at least one instance 
the arms are fully separated, VqT > R. 

(3) In the case of small interferometers the analytic ex- 
pression for the contrast degradation allowed us to sug- 
gest a simple recipe canceling the destructive effect of 
interactions completely. 

(4) In the case of large interferometers the effect of 
interactions can not, to our knowledge, be canceled en- 
tirely. Furthermore, the interatomic interactions set an 
universal limit on the number of atoms involved in the 
interferometric process: N <C KVq/qxd. Notice that this 
bound depends on the characteristics of the atom and 
the waveguide only (where the "beam-splitter velocity" 
Vq is supposed to be linked to the atomic transition fre- 
quency), while neither the timing of cycle nor the size of 
the atomic cloud enter. 

(5) Using the method developed we have analyzed the 
results of the recent JILA experiment on Michelson inter- 
ferometer on atom chip. In spite of comparable strength 
of interatomic interactions and stationary longitudinal 
trap present in the experiment, our results indicate a rel- 
atively weak role of interactions in fringe contrast degra- 



JILA experiment, N a = as shown ' 

(q\ theory, N & = as shown -\- 

theory, N a = 3.0xl0 3 

theory, HO trap alone 

theory, interactions alone, iV a = 3.0xl0~ 
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Propagation time, ms 



theory, N„ = 4.5x10 
theory, HO trap alon<j 
theory, interactions alone, N a = 4.5x10 




3 4 5 6 
Propagation time, ms 

FIG. 5: Fringe contrast ratio vs. duration of the interfero- 
metric cycle corresponding to the parameters of the JILA ex- 
periment with Michelson interferometer on a chip, with mag- 
netic gradient as the phase element, (a) Curves correspond 
to TV = 3. x 10 3 atoms, the as same 4 ms point in the experi- 
ment. The theoretical predictions for the actual experimental 
numbers of atoms at every run are also shown, (b) The same 
as (a), but for N = 4.5 x 10 4 atoms, where a reduction of the 
fringe contrast to 50% is expected. 



dation. This finding is can not be traced to any small 
parameter, but is a mere interplay of numerical prefac- 
tors. A moderate ten-fold increase (to N = 4.5 x 10 4 ) in 
the number of atoms will reduce, for large interferometer 
case, the fringe visibility to 50%, even without an addi- 
tional longitudinal trap present, reflecting the universal 
limit outlined above. 

(6) Michelson interferometer is a closed-loop white- 
light scheme by design; it is supposed to produce clear 
fringes even for input sources with a short coherence 
length. As one can see from the Fig.|2b), the mean-field 
pressure leads to two distinct effects. The first effect is 
the change in the relative momentum of the interferome- 
ter arms; this is what we addressed in the present work. 
The second effect is the distortion of the interferometric 
path, as a result of which the path becomes open, and 



7 



thus the interferometer becomes sensitive to the longitu- 
dinal coherence. For zero-temperature condensates this 
does not lead to any loss of contrast. At finite temper- 
ature the degradation due to the broken interferometer 
loop becomes relevant, and we are going to study this 
effect in the nearest future. 
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